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Sporulation in Bacillus subtilis serves as a paradigm for the development of two 
different cell types (mother cell and prespore) from a single cell. The mech- 
anism by which the two different developmental programs are initiated has 
been much studied but is not well understood. With the help of existing and 
new experimental results, a mathematical model has been developed that re- 
produces all published in vitro experiments and makes new predictions about 
the properties of the system in vivo. 
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Introduction 



Multicellular organisms achieve cell differentiation by selective expression of genes but the 
mechanisms that trigger such differential expression are difficult to define. Single-celled prokary- 
otes offer several experimental advantages for studies of cell differentiation; a good example is 
spore formation in Bacillus subtilis which is induced by nutritional stress. During sporulation 
the single cell organism differentiates into two cell types, a smaller forespore (prespore) and 
a larger mother cell, separated by a septum. The two different developmental programs are 
initiated by the transcription factor a F , which associates with the core RNA polymerase in the 
prespore but not in the mother cell. In the mother cell a F is trapped in a complex with SpoIIAB 
(hereafter referred to as AB). While the pathways that lead to o F dissociation are well-defined, 
it is not yet clear why a F association with the core RNA polymerase is restricted to the prespore 

m. 

Crystallography has demonstrated that one molecule of a F binds asymmetrically across the 
two subunits of the AB dimer 0. The binding of a F prevents access of SpoIIAA (hereafter 
referred to as A A) to one monomer of AB, but allows docking of A A on the other monomer 
0. Steric and electrostatic clashes are caused by the docking of AA, and a F is subsequently 
displaced from AB (HE). AA binds across the nucleotide-binding pocket of AB, and is phos- 
phorylated when this binding pocket contains ATP, yielding inactive phosphorylated AA (AA- 
P) and ABADP. ABADP can either be sequestered in an AA- AB-ADP complex or exchange 
ADP for ATP. Free AB (in both its ATP and ADP forms) can rebind a F , and sustained a F re- 
lease is thus expected to require reactivation of AA-P by SpoIIE-mediated dephosphorylation 
0. 

The restriction of a F release to the prespore has been ascribed to two aspects of asym- 
metric septation: i) the volume difference between prespore and mother cell, together with the 
accumulation of SpoIIE (hereafter referred to as HE) on the septum 0; ii) loss of AB in the 
prespore, caused by a transient genetic imbalance, as described below dEJ. Accumulation of 
IIE on (both sides of) the septum results in an increase in HE activity in the smaller prespore 
and reduced activity in the larger mother cell. It has been suggested that the increase in IIE 
activity may lead to sufficient dephosphorylation of AA-P to trigger a F release specifically in 
the forespore [6], and recent biochemical experiments indeed reveal that a F -AB complexes are 
highly sensitive to small changes in IIE activity 0. The resulting release of a F in the pres- 
pore might be reinforced by rapid degradation of unbound AB iflOl 151. which would prevent 
reassociation of a F and AB. Degradation of AB is not balanced by protein expression in the 
prespore because of a transient genetic imbalance; this arises from the fact that about 70% of 
the chromosome (including the spoIIA operon in which AB is coded) is initially excluded from 
the forespore and is imported into it only 10-20 minutes after septum closure ifTTI . A further 
proposed result of the genetic imbalance is the degradation and loss of a IIE inhibitor from the 
prespore, which would further accelerate a F release 0. The nature of such an inhibitor is, 
however, still elusive, and while the concentration of IIE on the septum and the degradation of 
AB in the forespore have both been shown to be important for successful sporulation, interfer- 
ence with either one of them does not abolish sporulation IHl[l2|. No model yet explains how 
these small and transient differences act together to maintain the release of a F for at least lh. 
Moreover, it is unclear whether these processes are sufficient for the initiation of sporulation or 
whether other, still unidentified, factors play a role Q- 

By employing a combination of experimental and mathematical modeling tools we derive 
the regulatory circuit that controls a F release. A large amount of information is available 
about the properties of this system, and most of the kinetic and thermodynamic parameters 
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are known. Our model quantitatively reproduces the system's known behaviour in vitro, in- 
cluding the sensitivity of AB -a F complexes to small changes in IIE activity and in AA con- 
centration. The model predicts that SpoIIE relocation to the septum is sufficient for prespore- 
specific formation of a F -RNA polymerase holoenzyme and that the required sensitivity stems 
from SpoIIAB allostery, which we confirm experimentally. Stochastic fluctuations in gene ex- 
pression are buffered by coupled translation of SpoIIAA and SpoIIAB combined with degra- 
dation of unbound SpoIIAB, while ATP consumption is limited by the formation of SpoIIAA- 
SpoIIABADP. The transient genetic imbalance plays no important role. 



Construction of the mathematical model 
The pathway that controls o F release 

The key steps and players involved in the regulation of a F can be summarized as follows. a F 
binds asymmetrically across the two subunits of the AB dimer [2 | and is displaced by steric 
and electrostatic clashes that are caused by the docking of SpoIIAA (AA) 101 El. AA binds 
across the nucleotide-binding pocket of AB, and is phosphorylated when this binding pocket 
contains ATP, yielding inactive phosphorylated AA (AA-P) and AB-ADP. AB-ADP can either 
be sequestered in an AA- AB-ADP complex or exchange ADP for ATP. Free AB (in both its ATP 
and ADP forms) can rebind a F , and sustained a F release is thus expected to require reactivation 
of AA-P by SpoIIE-mediated dephosphorylation Q. 

A model for the pathway that controls a F release 

During early attempts to develop a model of the system described in Figure 1 A in the main paper 
it became clear that some extensions would be necessary. One important extension, confirmed 
by new experiments is that AB, which is a dimer, binds AA with positive cooperativity (see be- 
low). This cooperative behaviour increases the number of reactions that need to be considered 
in the model. The shaded part of the scheme summarizes the interactions of one AB confomer 
and therefore includes only the well-known reactions that have been reviewed extensively (see 
Hilbert and Piggot (2004) and references therein). In brief, binding of a F to AB-ATP (61/62 
- all following numbers in round brackets refer to reactions in our model) leads to the forma- 
tion of a a F ■ AB-ATP complex lfT3ll . a F in this complex cannot bind to core RNA polymerase 
and is thus termed inactive lfT3lfT4ll . Subsequent interaction with AA leads to the formation of 
o F AB • ATP- AA (43/44), from which a F is rapidly released (79/80) O- Consistent with struc- 
tural information derived from crystallography Gl|4||, binding of two AA to the AB dimer while 
a F is still bound is impossible Q. The bound AA can be phosphorylated either immediately 
(185) or after a further AA has bound (189) fl4l . AA-P can be reativated by the phosphatase 
action of IIE (193-195; 0H6|). AB-ADP can either bind AA (9-12, 35/36,158) or rebind 
o F - either immediately (63/64) or after ADP- ATP exchange (140-143,170-173) Q. ADP- 
ATP exchange is also possible when AA is bound to AB-ADP (162,164) and dimers with no, 
or two different, nucleotides bound (19-26,53-56,69/70,86/87,92/93,156/157,178/179,186,191) 
also need to be considered. 

Given the allosteric behaviour of AB the scheme had to be extended to include also a second 
AB confomer (non-shaded part) as well as interconversions between AB states (94-137). This 
leads to a duplication of almost all the interactions mentioned so far. Only in the case of AA 
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and a F binding to nucleotide-free AB could the unfavoured conformation be neglected, since 
a comparison of the relevant parameters shows that these interactions would be of significance 
neither in the in vitro experiments nor in the sporulating cell. In the diagram, the conformation 
that binds AA with low affinity and ADP with higher affinity is depicted as circles (AB), and 

the other state as squares (AB). As is usual for descriptions of allosteric proteins (e.g. E3) 
these two states may be referred to as the tense (T) and the relaxed (R) state respectively. 

SpoIIAB is an allosteric protein 

Our experiments with fluorescence quenching, a Scatchard plot, and surface plasmon reso- 
nance (SPR) all provide evidence for the allosteric behaviour of SpoIIAB. The non-cooperative 
binding mechanism of published models failed to fit results from fluorescence quenching data 
obtained with AA and AB-F97W, a fluorescent derivative of AB CHI. When AB-F97W was 
mixed with AA, less than an equimolar AA concentration was needed to quench half the fluo- 
rescence (Fig. H]A). The effect was even stronger with the lower affinity AA mutant, AA-S58A, 
suggesting that binding of one AA to AB induces a conformational change which results in 
quenching of the fluorescence on the second non-bound monomer. A Scatchard analysis con- 
firmed the allosteric behaviour of AB (see Fig. 2A in the main paper). The experiments had to 
be carried out with the AA-S58A mutant to avoid AA phosphorylation by trace amounts of ATR 
(The AA/AB ratio is apparently overestimated by a factor of about 1.3, presumably because of 
errors in the determination of protein concentrations, and these probably explain the v > 2 
result). Moreover, SPR binding kinetics obtained with AB bound to the sensor chip are better 
fitted by an allosteric model (see below). Consistent with the hypothesis that AB is an allosteric 
protein, crystallographic structures of AB in complex with either a F or AA [121 51 have shown it 
to be a flexible molecule. The conformation of AB is sensitive to the type of ligand bound, with 
changes occurring within the nucleotide-binding region and in the angle between the subunits. 
However, our model assumes that both the a F - and AA-bound forms of AB are predominantly 
in the R-state, with the result that changes between the states would not be visible in the crystal 
structures because of low occupancy of the T state. 

Parameter determination 

Lack of quantitative data and experimental detail have often restricted the predictive power of 
mathematical models. This study has greatly benefited from the large amount of published 
experimental information and the availability of a well-established in vitro system. Values for 
essentially all the key parameters could thus be determined experimentally. The almost 200 rate 
constants included in the scheme were reduced to 30 independent rate constants by assuming 
that the kinetic constants of one AB monomer are not affected by ADP, ATP or AA being bound 
to the other monomer other than through a change in the relative concentration of the R and T 
states. 

The 27 independent kinetic constants for the in vitro model can be grouped into those char- 
acterising the AB-AA interaction (reactions 1-58), the AB-a F interaction in the absence (59- 
76) or presence (77-93) of AA, the interconversion between the R and the T state (94-137), the 
ADP-ATP exchange (138-182), and finally the phosphorylation (183-192) and dephosphoryla- 
tion (193-195) of AA. The values of many of these parameters have already been published. 
However, most of the on- and off -rates have been determined by SPR; this technique requires 
the proteins analysed to be attached to a sensor chip, which may alter their behavior. Fluo- 
rescence quenching experiments do not have this shortcoming, and a re-analysis of important 
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kinetic parameters was carried out by fluorescence spectroscopy with the help of fluorescent 
mutants [fl~8l IT5l . Constants that had been determined on the assumption that AB is not an 
allosteric protein were refined to fit the model. 

AB-AA binding constants 

The first group of constants are those characterising the AB-AA interaction (1-58). The rate 
constants for the formation of the AB-AA complex in the presence of ADP can be determined 
by SPR following the protocol established by Magnin et al (1997). These measurements can 
be carried out with either AB or AA bound to the sensor chip. However, AA-interactions 
with both AB monomers will occur only when AB is bound to the chip, since the spacing of 
chip-bound AA is unlikely to be appropriate for AB to bind to two molecules of AA. For the 
monomeric interaction when AA is chip-bound a total of 15 experiments were pooled from Shu 
et al (2004) together with other unpublished results from this laboratory. The on- and off-rates 
were determined as k* on = 1.6 x 10 5 M _1 s~ 1 (±2.1 x 10 5 standard deviation) and k* of j = 

6.85 x 10~ 3 s _1 (±4 x 10 -3 standard deviation) respectively. Here the off-rate will correspond 
to the the high affinity state, since upon binding of AA, AB will undergo a conformational 
change to the high affinity form. 

To facilitate the fitting of the binding of two AA molecules to one chip-bound AB molecule 
by means of the Biacore software we exploited the fact that one conformation of the AB and 
AB-AA 2 dimers is particularly favored; we thus simplified the binding reaction to a 3-step 
mechanism involving an initial binding of one AA to the AB dimer (in the low affinity state), 
a subsequent conformational switch and the binding of a second AA to AB-AA in the high 
affinity state. Fitting of SPR data with this 3-step mechanism was substantially improved over 
that with a single step (non-allosteric) model, with the average x 2 value of the fit being reduced 
by approximately 10-fold. (In contrast fitting was not improved by the 3-step mechanism when 
AA was bound to the chip.) The estimates obtained from 5 unpublished experiments for the on- 
and off -rates of the two binding steps were k on i — 1.1 x 10 6 M _1 s~ 1 (± 10 6 standard deviation), 
k on2 = 4 x 10 4 M-V 1 (±2.9 x 10 4 standard deviation), k ofn = 2.4 s" 1 (±2.2 standard 
deviation), = 8.4 x 10~ 3 s" 1 (±1.7 x 1CT 3 standard deviation). The conformational 

change occurs at 7.4 s _1 (±5.6 standard deviation) and 2 s _1 (±1 standard deviation) towards 
the high and low affinity states respectively. The affinity of the second bound AA will probably 
have been underestimated by SPR since attachment of AB to the grid is likely to obscure at least 
one binding site and to constrain conformational changes. Equally the AA-AB on-rate will have 
been underestimated when AA is fixed to the grid since AA, the smaller protein, will diffuse 
more rapidly when both are in solution. We therefore used an intermediate value k on i = 8 x 10 5 
M _1 s _1 as the general AA-AB on-rate in our simulations, and 7.4 x 10~ 3 s _1 as the off-rate 
in the high affinity state, which is the average value of all measurements made. and the 

rate for the conformational change in the disfavoured direction were used as determined by 
SPR above. In order to fit AA-AB binding as measured in fluorescence quenching experiments 
(Fig. [TJ the conformational change towards the favoured state had to be taken to be higher than 
estimated by SPR. We used 20 s" 1 ; a higher value would further improve the fit in Figure Hk. 
A likely explanation for the lower value for the SPR measurement is that fixation to the chip 
impairs AB conformational changes. 

The binding constants in the presence of ATP are more difficult to measure since AB-ATP-AA 
complexes are not stable because of phosphorylation of AA under these conditions. SPR mea- 
surements of the interaction of AB with the AA-S58A mutant of AA, which cannot be phos- 
phorylated by AB, revealed a higher affinity in the presence of ATP than of ADP [5 1. However, 
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the affinity of the AA-S58A mutant in the presence of ADP is lower than for the wildtype. 
We therefore could not directly use the values determined for complexes of AA-S58A with 
ABATP, and had to set the second off-rate to a lower, arbitrarily chosen value (k fj2 = 0.001 

In the absence of nucleotides AA-AB binding is weak Hl9\ and we used an off-rate that 
corresponds to a dissociation constant of 2.6 mM. 

It was assumed that the presence of a F did not affect the AB-AA affinity lfl~8l . 

a F -AB binding constants 

The next group of constants for which SPR data have been published Q are those that charac- 
terize the a F -AB interaction (59-93). While these could be used for the cr F -ABATP interaction, 
fluorescence spectroscopy data suggest a lower cr F ABADP affinity than that measured by SPR 
(Fig. EK)- Accordingly the a F ABADP off-rate was increased to 4 s _1 . To reproduce the kinet- 
ics of <j F -AB-ADP binding we have to assume that nucleotide-free AB can also interact weakly 
with a F (k ff = 40 s _1 ). The a F AB affinity in the presence of AA (77-93) was determined 
from fluorescence quenching data on AA-induced a F release. Fitting of the experimental data 
required the offrate with ATP present to be 100 times higher in the presence of AA than in its 
absence (a feature that had been noted previously lfT5l ). and required an even higher off -rate in 
the presence of ADP. We set the off -rate in the presence of ADP to the value of the on-rate, thus 
assuming no binding, but smaller values for the off -rate would give similar simulation results. 

Rates for AB conformational changes 

The SPR experiments on AB-AA binding provide the only rates available for the interconver- 
sion between R and T-states, and we used the rate for the disfavoured direction (2 s" 1 ) through- 
out (116-137). In order to obtain allosteric behaviour, the favoured direction (94-115) when 
both monomers have symmetrical ligands must be attained faster than the rate estimated by 
SPR in the presence of ADP when only one monomer has bound AA (7.4 s _1 ). We used 5 x 10 3 
s _1 . The exact value for this depends on the choice of the lid closure rate in the low affinity 
conformation (see below). 

ADP-ATP exchange 

The AB nucleotide binding pocket is partially covered by a flexible loop known as the "ATP- 
lid" 0. ADP-ATP exchange (138-182) requires lid opening lfl~8l and subsequent release of 
ADP and binding of ATP. The affinity of AB for ATP and ADP has been determined to be 
Kd ~ 200 /jM ll20ll . Fluorescence experiments suggest a slightly lower AB-ADP affinity 
lfT8ll . Assuming that the on-rate for this interaction is in the range found for other associations 
between proteins and small ligands (k on ~ 10 7 M _1 s _1 ifPTll ) this gives fc // ~ 2 x 10 3 s _1 . 
ADP-ATP exchange has been found experimentally to proceed at a rate of about 1 s _1 . This 
is much slower than would be expected from the nucleotide on- and off-rates, and has been 
interpreted as the lid opening rate [ 1 8 1 . The nucleotide-bound monomer is predominantly in the 
low-affinity conformation, and lid opening in this conformation was therefore taken to proceed 
at 1 s _1 . Given the choice for the rate of the conformational change, the lid closure rate had 
to be set to 5 x 10 5 s -1 in order to reproduce the observed rebinding rate of a F AB complexes 
after the AA-induced a F release (Fig. EtO, in agreement with the intuitive consideration that 
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the lid should be predominately closed when nucleotides are bound. It should be noted that 
these several steps associated with ADP-ATP exchange were included in the model only when 
the level of detail of available experimental information required a similar level of detail in the 
model - otherwise only the rate-limiting steps were considered in order to reduce the complexity 
of the model. 

AA phosphorylation and dephosphorylation rates 

The phosphorylation of AA (183-192) exhibits a biphasic timecourse (Fig. ElQ ||2T1 151. The 
initial rate is equal to the AA phosphorylation rate, and the subsequent reduction in rate is 
due to the sequestering of AA and AB in AB-ADP-AA complexes. The phosphorylation rate 
of AA can thus be estimated from the initial slope of the biphasic phosphorylation reaction 
as 0.013 s~\ a value that is close to previous estimates ll22l . The rate of AA-P production 
when AB-ADP-AA complexes have already formed is much higher than the rate of dissociation 
of the AB-ADP-AA complex. ADP-ATP exchange must therefore be possible while AA is 
bound to AB-ADP (162-165) - though at a reduced rate. Modelling of the steady-state rate of 
phosphorylation requires a ADP-ATP exchange rate for AB-AA of 1.2 x 10~ 3 s" 1 (Fig. EJC). 
ATP hydrolysis will be much slower than A A phosphorylation and was therefore ignored. 

Interestingly, when AB is pre-incubated for 5 minutes with AA and 5/iM ADP before ad- 
dition of 100 fjM ATP the biphasic phosphorylation timecourse can be directly transformed 
into a linear timecourse with the lower phosphorylation rate throughout (Fig. E]C); by contrast, 
pre-incubation with 5/iM ADP alone has only a minor effect (Fig. E]C). From these results it 
can be concluded that ADP-AB binding is weak but can be facilitated by AA, which itself binds 
weakly to AB when no nucleotide is bound |[T9l . In order to obtain such weak AB-ADP binding 
(while retaining a high affinity once ADP is bound), the R state must be favoured when the lid 
is open; in this conformation AB must have a lower affinity for ADP (k a ff = 10 5 s" 1 ) and the 
lid must be less mobile, that is the lid opening and closure rates need to be of order 10~ 4 s^ 1 . 
The only exception to this is when a F is bound, when we have to assume a lid closure rate of 
70 s -1 to reproduce a F -AB.ADP binding (Fig. EK)- In order to limit the number of variables 
we ignored the (disfavoured) closed lid state of AB when no nucleotide is bound; instead we 
assumed that the R state was 50 times more abundant than we had assumed elsewhere in the 
model. This is important in order to reduce the fraction of AB in the open, nucleotide-free 
tense state, which has a higher affinity for ADP than the relaxed state. Note that these particular 
parameter choices are needed to reproduce the in vitro results shown in Figure [2]but have little 
or no impact on the other results shown. 

Dephosphorylation of AA-P (193-195) in vitro with HE domain III has been measured to 
proceed at 0.085 s _1 ll23l . The on-rate for IIE-AA-P is unknown, as is the off-rate, and these 
were set to "standard" protein binding and unbinding rates, e.g. 10 6 M _1 s _1 and 0.1 s _1 respec- 
tively. 
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Verification of the model by comparison to in vitro data 

The model quantitatively reproduces experimental results of a F AB sensi- 
tivity to variations in AA and HE. 

With all the kinetic parameters determined we first sought to reproduce recent biochemical in 
vitro results (shown as circles in Fig. 0A) which revealed that a F ABATP complexes are highly 
sensitive to changes in HE activity [9|, thus lending support to the hypothesis that the septation- 
dependent increase in HE activity may lead to a F release in the prespore (6|. Translation of 
the scheme into a set of coupled differential equations allowed us to reproduce quantitatively 
the high sensitivity of <r F -AB-ATP complexes to such changes in HE activity, and also showed 
that the complexes are sensitive to changes in the AA concentration (compare 2.5 /jM and 4 
/iM curves in Fig. |3]A). The model also captured the kinetics of <r F -AB complex formation 
upon addition of ATP, the dissociation of the complex upon addition of AA and the effect of 
HE on re-formation of the complex (Fig. 2C in the main paper). The considerable variation in 
the experimental timecourses (compare the three black curves for IIE=0 nM in Fig. 2C in the 
main paper) may be due to the high sensitivity of the kinetics to concentration. Comparison to 
the results in Fig. |2j3 suggests that the slower recovery rate found in the simulation is the more 
realistic one. The values for the steady-state dissociation of the complex are more robust, and 
the standard deviations between three independent experiments, as indicated by error bars in 
Figure OK, are generally small. In both cases the modeling results are within the experimental 
range. The close fit of simulation and experiment, despite inherent experimental errors in the 
parameters including protein concentration, can be attributed to the robustness of the model (see 
below) and to the fact that the constants for the various on- and off -rates were all derived from 
similar experimental protocols. 

The robustness of the mathematical model to variations in parameters. 

The regulatory network can be regarded as robust if small changes in the parameter values 
lead to only small changes in the response. Given the availability of experimental data, the 
sensitivity of the model to variations in the kinetic parameters can be assessed by determining 
the deviation Aj of the simulation values s, from the TV experimental equilibrium results 
(Fig. |3K) in response to a change in the kinetic parameter value kj as 



Those reactions j, for which Aj did not increase by more than 1% when kj was either set 
to zero or varied between 1CT 3 and 10 6 in multiples of 10 were considered insensitive. Those 
for which even a small (10-fold) change in the parameter led to an increase in Aj by more than 
1% were considered sensitive. This analysis picks out nicely the key reaction pathways in the 
scheme and shows that relatively few reactions are sensitive to parameter variations. Sensi- 
tive reactions include the phosphorylation and dephosphorylation of AA, the binding of AA to 
a F AB complexes and also those that lead to ABADPAA complex formation or ADP-ATP 
exchange. 




(1) 
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Predictions of the effect of AB-ADP-AA complex formation. 

AB-ADP-AA complex formation coincides with a F release from the a F -AB-ATP complex 
(compare Fig. |3]A and B) and is important for the high sensitivity of a F release to changes 
in the AA concentration. The formation of AB-ADP-AA complexes also minimizes the need 
for cycling between AA and AA-P and associated ATP consumption. Thus, the ATP con- 
sumption per minute per a F released does not increase monotonically with increasing HE con- 
centrations, but reaches a minimum (Fig. |3] B, dotted line) when sufficient unbound AA is 
available to sequester AB in AB-ADP-AA complexes (Fig. |3]B, solid line). 

The AA-dependent trapping of AB-ADP in AB-ADP-AA complexes can be observed exper- 
imentally in the form of an increasing delay in a F -AB complex re-formation when equimolar 
concentrations of AB and a F are mixed with increasing molar ratios of AA (Fig. 2D in the main 
paper). The non-linear increase in the delay times with increasing AA concentrations reflects 
the non-linear AA-dependence of AB-ADP-AA complex formation. It should be noted that the 
extent to which this delay increases with increasing AA concentration would be smaller if the 
AA-AB-ADP affinity were lower. 



Application of the model to the sporulating cell in vivo 

With a model that reproduces the in vitro data quantitatively we are now in a position to in- 
vestigate the physiological situation. Relatively few alterations have to be made to the model 
to take account of the differences between the in vitro and in vivo conditions. The membrane 
location of HE in the cell does not need to be considered explicitly in the model since, given the 
dimensions of the forespore, spatial concentration gradients on the timescale of the reactions 
would be expected to be negligible. A change has, however, to be made to the ATP and protein 
concentrations, which are higher in vivo at the time of septation than in the in vitro experiments. 
The ATP concentration was set to 1 mM because the lower end of the experimentally reported 
ATP concentration range (0.85 mM to 3 mM) for vegetative B. subtilis |j24l|25l|26|) is likely to 
apply to the starving sporulating condition. The ADP concentration was set to 100 /iM, which 
is a tenth of the ATP concentration. To simulate the experimentally observed time course of 
protein expression during sporulation J5l|271|28l we started the sporulation reaction at time - 
120 minutes without protein and used an expression rate of 6 x M s _1 for AA and AB, 
whose translation has been suggested to be coupled [29), and 2 x 10~ 9 M s _1 for a F and HE 
(Fig. |4K). Septation (a 4-fold increase in HE and associated AA-P) was introduced after 120 
minutes 0|27l|28l unless otherwise stated. The transient genetic imbalance, which arises upon 
septum formation and which results in the prespore lacking the spoIIA operon for 10-20 minutes 
after septation 0IH1, was modelled by prohibiting expression of AA, AB and a F for 15 minutes 
after septation. The emergence of unphosphorylated AA (Fig. |4j3) further leads to repression 
of spoIIE and the spoIIA o pero n via inactivation of SpoOA IBOll . Note that in agreement with 
experimental observations 1311 almost half of all unbound AA (e.g. AA that is not bound to 
HE) is unphosphorylated already before septation (Fig. EJ3). 

An amino-acid sequence in the extreme C-terminus of AB markedly increases the AB degra- 
dation rate in vivo ifTOI . In the simulations unbound AB was therefore associated with a degra- 
dation rate of 4.1 x 10~ 4 s _1 , which gives rise to the experimentally observed half-life of about 
28 minutes ifTOl . 

The only rate constant that needs to be altered from the in vitro model is the HE phosphatase 
rate, which experiments suggest is lower in vivo ll32l . Modelling predicts that the reduction 
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needs to be by at least 10-fold to avoid septation-independent release of a F (Fig. 3 A in the main 
paper); a 20-fold reduction is used in all simulation results presented here. This reduction can 
be partly accounted for by the different in vivo and in vitro concentrations of Mn 2+ and Mg 2+ . 
The phosphatase rate is substantially higher in the presence of Mn 2+ ll33l : but, given that the 
phosphatase rate has been determined in the presence of Mn 2+ while Mg 2+ predominates in the 
bacterial cell, the rate in vivo may be considerably smaller than the literature value. However, 
this is unlikely to account for a 10-fold (let alone a 20-fold) reduction, and it may be that further 
alterations must be made to the HE environment or that an unknown HE inhibitor exists ll3Tll32l . 

HE accumulation on the septum, AB allostery and a short half- 
life enable compartment-specific gene transcription 

While septal accumulation of HE accumulation and a short half-life of unbound AB ifTOll 
have been established experimentally as key components of the mechanism that enables prespore- 
specific a F release, it is unclear how they contribute and whether they are sufficient given that 
the HE activity increase in the prespore is small and that the half-life of unbound AB strongly 
exceeds the time scale on which a F is activated. Moreover, it has so far been a paradox, how 
released a F can displace a A on the RNA polymerase despite the (unfavourable) competition of 
a F with the transcription factor a A (27). 

We will show in the following that septal accumulation of HE is sufficient because HE accu- 
mulates on the septum together with its substrate AA-P; this increases the rate of AA formation 
by at least 2.5-fold. The here-described allostery of AB enables the high sensitivity to such 
small changes in AA production; degradation of unbound AB is important to tune the protein 
concentrations in a way that small changes can act as a trigger without rendering the regula- 
tory network vulnerable to fluctuations in protein expression. The relative affinities of the RNA 
polymerase for a F and a A are shown to be optimized to allow for a coexistence of both holoen- 
zymes. Other hypotheses such as the removal of a HE inhibitor or effects of the transient genetic 
imbalance are shown to be irrelevant. 

Accumulation of HE (in association with its substrate AA-P) at the septum 
is sufficient to induce o F release in the prespore. 

Contrary to previous suggestions, the model demonstrates that septal accumulation of HE is 
sufficient to induce o F release in the prespore; it is not necessary to increase the HE phosphatase 
rate in the prespore (for example by the removal of an inhibitor). As AA production is equal to 
the phosphatase rate times the concentration of HE- AA-P complexes, it is sufficient to increase 
the substrate and enzyme concentrations by accumulation on the septum (Fig. EK)- Substrate 
and enzyme increase together when HE accumulates on the septum because, owing to low 
phosphatase rate and the slightly higher total AA concentration relative to the HE concentration 
in the cell, almost all HE is associated with AA-P before septation (Fig. 3B in the main paper); 
enzyme and substrate would accumulate together on the septum even if the IIE-AA-P affinity 
were as low as micromolar. 

An increase in either the HE activity or AA-P concentration alone is far less efficient even if 
the phosphatase rate is set to the highest physiological possible value (Fig. |5] A). It is therefore 
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necessary for HE to be present before septation, since that is the only way in which septal ac- 
cumulation of IIE could also increase the AA concentration in the prespore. In fact asy mmetric 
septation appears to require the presence of a significant IIE concentration Il34ll35ll36ll37ll. 

Formation of the cr F -RNA polymerase holoenzyme 

Prespore-specific a F release is only a first step to differential gene expression. In order to direct 
gene expression a F needs to bind to the RNA polymerase, from where it needs to displace a A . 
The mechanism of the latter has so far been unclear given the (unfavourable) competition of 
a F with o A {a A binds to the polymerase with a 25-fold higher affinity (560 nM and 22.5 nM 
for a F and a A respectively), but it is only slightly less abundant (~ 7.5// M) than released 
a F E3). As the core RNA polymerase concentration is comparatively high(~ 7.5// M), the 
mathematical model predicts that micromolar concentrations of a F -RNApol holoenzyme can 
form spontaneously in response to asymmetric septation (Fig. 3E in the main paper), and that 
these complexes co-exist with micromolar concentrations of cr^-RNApol holoenzyme (Fig. 3E 
in the main paper). The latter result is in agreement with the experimental observation that 
a^-RNApol holoenzyme persists throughout the sporulation process ll38l . 

Figure 3E in the main paper further shows that IIE (and AA-P bound to it) needs to increase 
by at least 2.5-fold in order to ensure the formation of micromolar cr F -RNApol concentrations. 
Given the difference in volume between the two compartments, accumulation of IIE on the 
septum will in fact lead to at least such an increase (see Supplementary Information). However, 
as discussed in the Supplementary Information, a 4-fold increase is a more realistic description 
of the sporulating cell and is therefore used in all the simulations presented here. 

The role of the transient genetic imbalance 

The view that a F release is driven by a change in IIE activity at the asymmetric septum has 
been challenged by the finding that a F was still released, albeit with a significant delay, in 
mutants where wild-type IIE was replaced by a soluble phosphatase domain (IIEAmem) II 1 21 . 
The formation of cr F -RNApol in the IIEAmem mutant was subsequently shown to be enabled 
by the block in AB expression that follows septation due to a transient genetic imbalance 

The simulation correctly predicts the phenotype of mutants that lack the transient genetic 
imbalance and the impact of the IIEAmem on such background (Fig. EDA)- Importantly however, 
the simulation also shows that a 10-20 minute block in AB expression (as may arise from 
the transient genetic imbalance - Fig. |5jB) is insufficient to drive prespore specific a F release 
in IIEAmem mutants (Fig. compare solid and dotted lines). A continuous block in AB 
expression is necessary (Fig. |6£, solid line) and this cannot be achieved by SpoOA repression 
since, unlike wildtype, IIEAmem does not dephosphorylate sufficient AA upon septation (Fig. 
03). 

We conclude that a genetic imbalance as short as 10-20 minutes cannot drive compartment 
specific gene expression. AB decay leads to o F release only after 30-40 minutes and therefore 
only enables significantly delayed a F release in IIEAmem mutants. Successful sporulation, 
however, requires that a F is released within 10 minutes after septation [1J. IIE accumulation on 
the septum is therefore necessary to enable rapid a F release. Additional mechanisms such as 
inhibitor removal or a block in the release of dephosphorylated AA from IIE before septum 
formation iTTTll are not required. 
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The role of the short half-life for unbound AB 



While results like those above question a role for a short half-life of unbound AB in rapid 
prespore- specific a F release, mutants with stable AB are sporulation deficient ifTOl . The sim- 
ulation now reveals that the 28 minutes half-life of unbound AB is important to limit the AB 
concentration before septation (Fig. 3D in the main paper). The use of AB degradation rather 
than a lower expression rate is important to guarentee robustness to stochastic fluctuations in 
gene expression as is discussed below. 

Allosteric behaviour of AB ensures both sensitivity and robustness 

AB allostery is key to the high sensitive to changes in the IIEAA-P concentration and ensures 
that the system is robust to the initial parallel increases in the protein concentrations (Fig. |7K)- 
Due to AB allostery <r F -AB complexes are stable at low AA concentration but become rapidly 
dissolved once the AA concentration reaches a critical threshold after septation. The premature 
formation of AB-AA complexes in the absence of allostery (Fig. |7j3) would have two detri- 
mental effects. First, sequestration of AB would free a F (Fig. |7K,B)- Secondly, sequestration 
of AA would reduce the fraction of IIE that is bound by its substrate AA-P (Fig. |7JZ!), thereby 
interfering with the simultaneous accumulation of IIE and AA-P when the septum forms, and 
thus limiting the formation of unphosphorylated AA in the prespore since - unlike in metabolic 
pathways - IIE's substrate (AA-P) is a protein and present at a low concentration in the cell. 

We have investigated whether non-allosteric behaviour of AB could be compensated for by 
a stronger inhibition of IIE. However, the model shows that even a 1000-fold inhibition of IIE 
activity would be ineffective in preventing premature a F release (Fig|7J dashed lines) and would 
fail to provide for specific a F release in response to septation. 

Economic efficiency through limiting ATP consumption 

The avoidance of "ATP wastage" can be expected to become particularly important once nutri- 
ents are scarce as is the case at the onset of sporulation. The regulatory network that controls a F 
release is optimized to limit ATP wastage. Rapid formation of a F -RNApol holoenzyme (Fig. 
0\) is accompanied by the accumulation of most of the AB in AB-ADP-AA complexes (Fig. 
|5jB). If the simulation is altered so as to prohibit the formation of these complexes, AB accumu- 
lates instead in ABATPAA complexes; but the consumption of ATP is then greatly increased 
(Fig. 4D in the main paper). We conclude that, for a F -RNApol holoenzyme to be formed and 
to persist without excessive waste of ATP, the formation of AB-ADP-AA complexes and the 
degradation of AB are both important. AB-ADP-AA complexes therefore play an important 
role in acting as a sink for AB after the increase in IIE (and associated AA-P) that occurs at 
asymmetric septation. 



Robustness 

Biological networks need to be organized in a way that they are sensitive to regulatory input but 
robust to the small variations that may arise from changes in environmental conditions or from 
stochastic fluctuations within the cell. The sporulation network meets both conditions to an 
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impressive degree. The protein concentrations are tuned to ensure a high sensitivity to changes 
in the IIE concentration while being robust to stochastic fluctuations in gene expression rates 
(Fig. 4A in the main paper). The additional robustness of the network to random changes 
in values in kinetic parameter and enzymatic rates (Fig. |3) justifies our approach of using 
the kinetic parameters that were determined from in vitro experiments. By means of control 
coefficients we determine the reactions that the regulatory network is most sensitive to as the 
phosphorylation and phosphatase steps. 

Protein expression is tuned to ensure robustness 

During the vegetative phase, the concentrations of a F , AA, AB, and IIE are too low to allow 
sufficient a F -RNApol formation to drive prespore development (Fig. |4j\ and Fig. 4B in the 
main paper). This efficiently protects the cell from premature a F activation. The triggering 
of differential gene expression requires the formation of IIE-AA-P complexes, and therefore 
the expression of AA and IIE, prior to septation (see above). In order to avoid premature a F 
activation all sporulation proteins (a F , AA, AB and IIE) need to be expressed in parallel prior 
to septation. Given the low numbers of proteins expressed, the process is inherently stochastic. 
The ratio between the proteins is essential to avoid spontaneous a F release, and translational 
coupling of all proteins is the best strategy to avoid failure (dashed line in Fig. 4A in the 
main paper). However, the downside of such coupling is that the ratios between proteins would 
then need to be established by degradation which is uneconomic. The simulation reveals that 
already the coupling of AA and AB translation strongly increases robustness as compared to a 
mechanism that does not involve coupling (compare dashed and dotted lines in Fig. 4A in the 
main paper). 

The robustness will even be higher than suggested by Figure 4A in the main paper if the 
time point of septation depended on the protein (and in particular the IIE) concentration as may 
indeed apply to the cell Il35ll37ll . The impact of the timepoint on sporulation efficiency reflects 
the need to exceed a certain protein concentration before the septation-dependent increase in 
IIE-AA-P can trigger the formation of a F -RNApol at micromolar concentrations. The dotted 
line in Figure 4B in the main paper suggests this threshold to be 10 fi M for a F (and corre- 
sponding concentrations for the other proteins) at the timepoint of septation which is within the 
physiological concentration range for which a F -RNApol formation has been observed II51 I281 . 
We suggest that in order to avoid unproductive septation, the timing of septum formation must 
be linked to the overall concentration of the proteins relevant to sporulation; this may explain 
the large variance in the delay between the onset of sporulation and septation that is observed 
under different sporulation conditions. In the absence of a septation-dependent increase in 
the IIE-AA-P concentration the concentration of o" F -RNApol remains too low to trigger down- 
stream signaling (solid line in Fig. 4B in the main paper). This protects the mother cell, where 
the concentration of IIE-AA-P actually decreases slightly, from spontaneous a F release. 

Critical parameters that determine the behaviour of the regulatory net- 
work 

The impact of the individual kinetic parameters on the behaviour of the model can be assessed 
by means of control coefficients, as previously suggested for the analysis of signaling networks 
by Lee and co-workers ll39l . Control coefficients can be used to quantify the relative change in 
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a model variable in response to (small) relative changes in a parameter value h. Here they were 
determined from 

r RNA P oi(a F ) _ SRNApol(a F ) kj 
* ~ RNApol(a F ) Ski' 

allowing the impact of parameter perturbations on the concentration of RNApol(a F ) com- 
plexes to be ascertained. To understand the relative impact of the different interactions we 
studied the 30 independent rate constants listed in Tables EKd Control coefficients were de- 
termined in the equilibrium state before septation. We used 15 fjM <r F , 30 fjM AA, 22.5 /xM 
AB, and 15 /iM HE for the concentrations of the sporulation proteins, which reflect the exper- 
imentally observed ratios between them and are within the experimental concentration range 
10 |27l |28l. As required by the summation theorem for concentration control PRl they sum 

to zero. Several reactions exert strong control ^(j^ NA P° l( - a ) > 1^ m0 st notably AA phospho- 
rylation and dephosphorylation, but the reactions responsible for ADP-ATP exchange and the 
binding of the transcription factors also have a strong effect (Tables [T^ISl*- 

Affinities of the transcription factor are particularly interesting since the relative affinities of 
a F and a A appeared counterintuitive when first determined ll27ll . The mathematical model now 
reveals that an equal, or higher, affinity of a F relative to a A for the RNA polymerase would 
result in a strong increase in RNApol(cr F ) complex formation even before septation, which 
could lead to premature activation (Fig. 3D in the main paper). 

We conclude that several safeguards are employed to avoid premature activation of sporu- 
lation notwithstanding the high sensitivity of the regulatory system to small perturbations in 
protein concentration. These safeguards include a tight control of protein concentrations, an 
adjustment of the phosphatase rate, the different affinities of the RNA polymerase for the two 
transcription factors and the allosteric behaviour of AB. 



Discussion 

Complex biological networks are difficult to understand from verbal descriptions, and mathe- 
matical modelling is a promising tool to integrate the available information. However, realistic 
models cannot be built without quantitative experimental data. The system studied here offered 
several unusual advantages in that the entire regulatory network could be reconstituted in vitro 
and specific fluorescence probes were available to study all relevant interactions. This enabled 
realistic values for all kinetic parameters that affect the model to be estimated from experiments, 
and allowed the derivation of a very detailed model of the process that reproduced the experi- 
mental data quantitatively. In applying the model to the physiological situation the phosphatase 
rate of HE, which exhibits a lower rate in vivo than in vitro B32I . was the only parameter that 
had to be changed. 

We began with a simple, basic circuit from published schemes but this rapidly grew into a 
large network after mathematical modeling predicted that AB is an allosteric protein. Here the 
combination of theory and experiment became particularly fruitful because theoretical predic- 
tions led to new experiments and analysis. Some rate constants, such as the phosphorylation rate 
and the <r F -AB-ADP binding rates, which had been estimated from previous experiments using 
simple models, also needed to be re-evaluated as the regulatory scheme became more complex. 
Once the rate and binding constants had been determined, the time-dependent variation in the 
system variables (e.g. concentrations of the pathway components and their complexes) could 
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be analysed by translating the interactions into a set of ordinary differential equations, which 
could be solved numerically. Comparison of the model predictions to data derived in vitro con- 
firmed that the model was realistic and a perturbation study revealed that the model is robust to 
variations in the kinetic parameter values which justified a direct application of the model to the 
physiological situation. 

Extrapolation of the model to the physiological situation enabled us to evaluate many effects 
that had been associated with the regulation of a F release in the prespore and to integrate them 
into a mechanistic picture of the regulatory process. We conclude from analysis of the model 
that HE accumulation on the asymmetrically placed septum is sufficient to trigger prespore- 
specific o F release because IIE association with its substrate AA-P maximises AA production 
in the prespore while AB allostery generates the required sensitivity. Thus, the volume differ- 
ence between the two compartments (mother cell and prespore), AB allostery together with the 
architecture of the regulatory network as well as the balanced concentrations of all sporulation 
proteins are the key factors that ensure compartmentalized gene expression during sporulation 
in B. subtilis. Alternative mechanisms such as a block in AB expression or the removal of 
an (elusive) IIE inhibitor, possibly through a transient genetic imbalance, are unlikely to be of 
major physiological relevance. 

The regulatory system that leads to compartmentalized gene expression in B. subtilis is 
optimized in many ways. First of all the effect of the septation-induced difference in the ratio 
of surface area to volume is maximized between the two compartments by concentrating the 
phosphatase IIE on the asymmetric septum, that is on the common interface. Given the small 
size of the bacterium, a homogeneous distribution of IIE on the surfaces of both compartments 
would have led to only a small change in the relative phosphatase activities. Association of IIE 
and AA-P ensures a simultaneous increase in concentration of both enzyme and substrate in the 
prespore, which in turn ensures the strong increase in AA production that leads to a F release. 
Secondly, the protein concentrations are tuned such that minute changes in IIE and AA can tip 
the balance while robustness to fluctuations in gene expression is ensured by coupled expression 
combined with protein degradation. Thirdly, the previously unnoticed allosteric behaviour of 
the AB dimer enables the high sensitivity to changes in IIE and AA and at the same time 
prevents premature AB-ADP-AA complex formation, which would lead to o F release. Finally, 
the formation of AB-ADP-AA complexes after septation serves as sophisticated mechanism 
for preventing a F from rebinding to the released AB, thus providing rapid and prolonged a F 
release while minimizing cycling and associated wastage of ATP. AB-ADP-AA complexes thus 
serve to remove AB (AB sink ||4ll|42|) rather than AA (AA sink ll43l0 . 

The model's prediction that pre-septational AB-AA complex formation is detrimental and 
that AB-ADP-AA complexe serve as AA sink stands in direct opposition to the suggestion that 
AB-ADP-AA complexes may act as an AA sink to prevent premature a F release ll43ll . This 
suggestion had been proposed to explain the counterintuitive experimental observation that un- 
phosphorylated AA can be observed long before septation but without inducing <x F -dependent 
gene transcription ll3Tll . Our model explains the pre-septational accumulation of (inactive) AA 
without invoking AB-ADP-AA complex formation (Fig. |4jB,|33). The large fraction of A A that 
is unphosphorylated at the very onset of sporulation (Fig|4j3) is the consequence of the low 
protein concentrations which prevent the formation of protein complexes; lack of AB-ATP-AA 
complex prevents AA phosphorylation. Similarly, most of the a F is free at that time, but its 
concentration is too low to allow the formation of sufficient <r F -RNA polymerase holoenzyme 
to effect gene transcription. Once protein concentrations increase, the model predicts that much 
AA will be bound in complexes with AB-ATP. Experiments find that half of AA is unphos- 
phorylated aready before septation, and the simulations predict that this applies indeed to the 
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soluble AA pool (not HE bound). 

The critical importance of the simultaneous accumulation of IIE and its substrate AA-P as 
well as the need to render septation dependent on the protein (and in particular IIE) concen- 
tration would probably have remained unnoticed without mathematical modeling. Much the 
contrary it had already been reasoned that the release of the product AA from IIE was inhibited 
before septation ll3T1l . 

The mathematical model also resolves the paradox of why degradation of unbound AB is 
necessary for sporulation while the half-life (~ 28 minutes) strongly exceeds the time scale on 
which o F needs to be released (~ 10 minutes) ifTOl [Til. The model reveals that the observed 
half-life reduces the AB concentration, whose expression is presumably coupled with AA ll29l . 
to a level that small changes in the AA concentration can trigger a F release. Coupled expres- 
sion and degradation of unbound AB (rather than all AB) confer robustness to fluctuations in 
gene expression. AB allostery ensures robustness to parallel changes in protein concentration 
(without this feature the initial increase in protein concentration might cause a F release in the 
entire sporangium). A IIE dependency of septation is likely to ensure that sufficient protein 
has been expressed such that septation can trigger the formation of micromolar a F -RNApol 
holoenzyme concentrations. The relative affinities of the transcription factors a F and a A for the 
core RNA polymerase are tuned to prevent premature activation whilst enabling the system to 
be triggered by small perturbations. Moreover RNApol(cr A ) holoenzyme remains present when 
RNApol(a F ) emerges, in such a way that general transcription programs can still be carried 
out. The robustness of the regulatory circuit to alterations in concentration and values of kinetic 
parameters is such that the process can be deregulated only by mutation. 

We conclude that the regulatory circuit described here appears to be optimized to the needs 
of the organism. It is tuned to ensure both high sensitivity to stimuli and robustness to random 
fluctuations. We show that the dimeric, allosteric kinase SpoIIAB plays a crucial role by en- 
abling sensitivity, robustness, and minimal ATP consumption. The high sensitivity to changes 
in the IIEAA-P concentration enables the two different developmental programs to be trig- 
gered by small concentration changes (such as can be induced by asymmetric septation) and 
with minimal ATP consumption. 

While asymmetric cell division, and the resulting difference between the two compartments, 
is an elegant mechanism for achieving compartmentalized gene expression, other symmetry- 
breaking mechanisms must exist to induce compartmentalized gene expression in those sporu- 
lating bacteria, such as Sporosarcina ureae, that divide symmetrically during sporulation P4l . 
A deeper insight into such alternative mechanisms will be of great interest. 
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Appendix: Comments on how to estimate of septation-dependent 
HE activity change 

Given that IIE accumulates on the septum (together with any AA-P that is bound to it) the 
number of IIE-AA-P complexes acting per unit volume increases in the prespore and decreases 
in the mother cell. A lower limit for the septation-dependent increase in IIE activity in the 
prespore can be determined by assuming a) that all the IIE locates, together with its substrate, 
on to the septum in such a way that its activity is displayed equally on both faces and b) that 
the cell is cuboid.. Insertion of a septum at an xth of the length of the cuboid will give rise to 
two compartments, with the smaller having an xth of the total volume. If half of all IIE-AA-P 
is present in the smaller volume then there is an effective |-fold increase in IIE-AA-P per unit 
volume of the smaller compartment. Septa form at one-fifth or less of the length of the cell P31l 
and the minimal increase in IIE activity is therefore 2.5. 

In fact, each of the assumptions mentioned above underestimates the probable increase in 
IIE activity: the total quantity of IIE is slightly greater on the prespore face than on the mother- 
cell face of the septum ll46l 1311 : and the prespore is not a cuboid but has rounded ends H71 . 
Thus, although a 2.5-fold increase in IIE activity is sufficient to ensure the formation of enough 
a F -PvNApol to permit transcription in the prespore (Fig. |5JC), the 4-fold increase that we have 
assumed throughout this paper is a more realistic representation of the physiological situation. 
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Figures 




Figure 1: SpoIIAB is an allosteric protein. (A) Fraction of AB (1 /xM total) unbound as a function 
of the concentration of AA in the presence of 100 /jM ADP, as determined by fluorescence quenching 
experiments with fluorescent mutant AB-F97W lfl~8l . The points show previously unpublished exper- 
imental results; the line represents the binding curve expected from the kinetic constants employed in 
the mathematical model (which includes allostery of AB). (B) Scatchard plot obtained by determining 
the fraction of bound AB (u) when 4/iM AA-S58A is mixed with different concentrations of AB. The 
points show previously unpublished experimental results; the lines represent the predicted Scatchard plot 
based on the assumption that AB is an allosteric protein (continuous line) or non-allosteric (dotted) with 
an affinity for AA-S58A-AB 3.7 times less than for AB-AA (Table QJ. Note that a 1.3-fold higher AB 
concentration than measured had to be used to keep v < 2. 
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Figure 2: Estimate of parameter values. (A) Interactions between 1.3 /jM g f and 1 fjM AB. The 

results of fluorescence quenching experiments [9] are shown as output from a recorder, and simulations 
are shown as smooth lines, either dashed (for ATP-dependent interactions) or continuous (for ADP- 
dependent interactions). (B) a F release from ABATP upon addition of AA. The results from two exper- 
iments are shown for the disruption of a F -AB complexes by 0.4 (red); 0.7 (yellow); 1 (blue); 1.5 (cyan); 
2 (magenta); 3 (black) AA. The corresponding simulations are shown in smooth black lines, which 
are labeled with the relevant concentrations of AA. Experiments were carried out as described in [ 9 1 with 
the fluorescent mutant cr F -W46L. (C) Phosphorylation of AA by AB when AB is directly incubated with 
40 fjM AA and 100 fjM ATP (o) or pre-incubated for 5 minutes with either 5 ADP (•) or 5 ADP 
and 40 fjM AA (□). The symbols show the experimental results fTJ], normalised by the factor 1.3 to 
take account of inaccuracies in the protein determination (see main text). The lines show the predictions 
of the mathematical model. 
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Figure 3: The dissociation of a F AB complexes is sensitive to changes in the concentration of 

HE and AA. (A) Steady state after IIE-dependent dissociation of ct f ABATP (1.3/iM a F , 1/xM AB 
dimer) in the presence of 2.5/iM (blue) or 4/jM (red) AA. Circles represent experimental results (9j- 
Continuous lines show the predictions of the mathematical model, which includes the formation of the 
AB ADP AA complex. Dashed lines show simulations based on the assumption that the formation of 
such a complex is impossible. (B) Kinetics of IIE-dependent dissociation of u^-AB-ATP (1.3/iM a F , 
1/jM AB dimer) in the presence of 2.5 /jM AA, in the absence of HE (black lines) or with 10 (green), 
40 (red) or 100 (blue) nM HE. Experimental results are taken from [9|, and are shown as output from 
a recorder. Note that three experiments were run in the absence of HE. The smooth continuous lines 
show simulations from the mathematical model. (C) Simulations of the proportion of AB sequestered 
into ABADPAA complexes (continuous line) and ATP consumption per minute per a F released (dotted 
line) in the presence of 2.5 /jM (blue) or 4/iM (red) AA. (D) er^-AB complex formation on mixing of 1 
fjM AB and 1.3 /xM a F with no (black), 2 fjM (blue), 4 /iM (green), or 6 fjM (red) AA. Experimental 
results (two for each concentration of AA) are taken from [ 9 1 and are shown as output from a recorder; 
the dashed lines show simulations from the mathematical model at the same concentrations of AA. 
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Figure 4: The simulation reproduces the experimentally observed concentrations of sporulation 
proteins. (A) Increase in the concentration of sporulation proteins during sporulation (AA- dashed; AB - 
continuous black; HE - dotted; a F - grey). (B) Time-dependent fraction of AA that is unphosphorylated 
(solid line), unbound and unphosphorylated (dotted line) or bound in HE AA-P complexes (dashed line). 
In both panels, sporulation-dependent protein expression is started at -120 minutes; a four-fold increase 
in HE and associated AA-P (to simulate the effect of asymmetric septation) is introduced at t=0. (C) 
Fraction of a F released in simulations based on the assumption that the HE activity in vivo was not 
inhibited (continuous grey line), or inhibited by 10 (dotted), 20 (continuous black) or 40 times (dashes) 
compared with the maximum HE activity measured in vitro. A four-fold increase in HE and associated 
AA-P (to simulate the effect of asymmetric septation) is introduced at t=0. 
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Figure 5: Simulations showing that a septation-dependent increase in HE AA-P is sufficient for 
a F release. (A) a F release in response to a 4-fold increase in HE alone (dashes) in AA-P alone (dotted), 
or in HE and AA-P (continuous line) imposed at t=0. The HE phosphatase rate was either 10-times (grey) 
or 20-times (black) smaller than in vitro. (B) Time-dependent fraction of HE in IIE-AA-P complexes. 
(C) Time-dependent a F release in wildtype (black), when AB expression is reduced by 30% (red line), if 
HE expression is increased by two-fold (green and blue lines) and if HE activity is increased by two-fold 
(green line). In all panels, the septum was assumed to form at t=0 (see main text). 
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Figure 6: The effect of a transient genetic imbalance on a F release. (A) Formation of RNApol(c F ) 
complexes, (B) total AB concentration for wt (continuous line) spoIIAB(ori) (dashes) IIEAmem & spol- 
/AB(ori) (dotted). (C) Formation of RNApol(cr F ) complexes, (D) total AA concentration for wt (con- 
tinuous grey line), for HEAmem if chromomosome translocation is blocked (solid line) if SpoOA is 
constituively repressed (dashed line) or if IIEAmem & spoIIAB(on) (dotted). In all panels, the septum 
was assumed to form at t=0 (see main text). 
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Figure 7: The allosteric behaviour of AB is necessary for successful sporulation. (A) The formation 
of RNApol(cr F ) holoenzyme, (B) the fraction of AB in AB ADP AA complexes, (C) the fraction of HE 
in IIE-AA-P complexes, and (D) the fraction of unphosphorylated AA if AB is allosteric (continuous 
line) or not allosteric and the activity of HE is inhibited 20-fold (dotted) or 1000-fold (dashes) compared 
with the maximum activity measured in vitro. A four-fold increase in HE, together with its associated 
AA-P, is imposed at t= 0. 
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Figure 8: AB ADP AA complexes form upon asymmetric septation. The simulations show the 
effect of a four-fold increase at t=0 min of HE, together with its associated AA-P, on (A) formation of 
RNApol((T F ) holoenzyme, (B) fraction of AB in AA2-AB ADP complexes. (C) ATP consumption per 
minute per a F released for wildtype (continuous line) or for a variant in which AB ADP AA cannot be 
formed (dotted). 
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Figure 9: The robustness of the regulatory network. ) (A) The fraction of successful sporulation 
events if the rates of protein expression are varied randomly by the factor given on the x-axis and all 
expression rates are coupled (-), AA and AB expression is coupled (•■■)> a U expression rates are varied 
randomly, the AB expression rate is 4 x 10 -9 M s _1 and there is no AB degradation (- -) . A success- 
ful sporulation event was defined as one for which the concentration of RNA polymerase bound to F is 
lower than 0.4 before, and higher than 1 fjM after, septation. Standard deviations are based on 10 
sets of 100 runs each. (B) Fraction of successful sporulation events if either all binding rates (dashed), 
all enzymatic rates (dotted), all binding and enzymatic rates (solid) are varied randomly by the factor 
given on the x-axis. Standard deviations are based on 10 sets of 100 runs each; successful sporulation 
events are defined as those with a concentration of RNA polymerase(<r F ) holoenzyme less than 0.4 fjM 
without, and greater than 1 fiM with, a 4-fold increase in HE (and associated AA-P). (C) Concentration 
of RNA polymerase(<r i? ) holoenzyme as a function of the reference protein concentrations (as exempli- 
fied by the a F concentration). The relative concentration of HE (and the associated AA-P) was either 
increased by four-fold (continuous line) or left unchanged (dotted). (D) Concentration of RNA poly- 
merase a F holoenzyme if the affinities of the core RNA polymerase for a F and a A were as in the wild 
type (continuous line) or made equal by varying the rate of dissociation of either the a F holoenzyme 
(k ff{a F ) = 0.02 s" 1 - dashes) or the o A holoenzyme (k ff(a A ) = 0.55 s" 1 - dotted). 
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Tables 



Table 1 : Control coefficients for the kinetic rate constants employed in the in vitro and 
in vivo simulations. The determination of control coefficients Q^ NA P° l ( a > [ s described in the 
main text. References to a particular figure show that a parameter value was derived from results 
presented in that figure rather than from previously published results. 



rates 


f~iRN Apol(cr' ' ) 


AA-AB on-rate 


2 x 1(T 5 


AA-AB off -rate 


-0.07 


AA-AB. ATP off-rate 


0.03 


AA-AB. ADP off-rate 


-0.02 


AA-AB off-rate 


5 x 10~ 7 


a* -AB on-rate 


-0.66 


cr^-AB ATP off -rate 


0.29 


o^'-AB ADP off-rate 


10" 3 


o-' -AB- off-rate 


-3 x 10~ 4 


a" -AB.ATPAA off-rate 


0.15 


a" -AB.ADPAA off-rate 


10" B 


AB conformational change 


-0.84 


AB conformational change 


0.08 


AB conformational change 


-0.08 


lid opening T state 


-1.16 


lid opening R state 


-9 x 10" B 


lid closure T state 


1.14 


lid closure R state 


4 x 10~ 7 


lid closure R state (a 1, bound) 


3 x 10~ 4 


nucleotide-AB on-rate 


-5 x 10~ 4 


nucleotide-AB off-rate 


-0.32 


ADP-AB off -rate (T state) 


-3 x 10" a 


ATP- ADP exchange for AB AA 


-6 x 10" 3 


AA phosphorylation 


-1.79 


AA-PIIE on-rate 


0.035 


AA-PIIE off-rate 


-0.034 


AA dephosphorylation 


see Table [2] 
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Table 2: Kinetic rate constants only employed in the in vivo simulation. The determination 



of control coefficients C- po is described in the main text. 



rates 


^RNApolia 1 ") 
°7 


simulation value 


Ref 


AA dephosphorylation 


3.24 


3.2 x 10-* 


Fig.|3\, [32| 


(7-RNApol on-rate 


0.33 


10 B 


|27| 


a" -RNApol off -rate 


-0.78 


0.55 


|27| 


cr A -RNApol off -rate 


0.45 


0.02 


|27| 


AB degradation 


n/a 


4.1 x 10~ 4 


HOI 



34 



